Enhanced resolution profiling in twins reveals differential methylation signatures of type 2 diabetes with links to its complications.

BACKGROUND: Type 2 diabetes (T2D) susceptibility is influenced by genetic and environmental factors. Previous findings suggest DNA methylation as a potential mechanism in T2D pathogenesis and progression. METHODS: We profiled DNA methylation in 248 blood samples from participants of European ancestry from 7 twin cohorts using a methylation sequencing platform targeting regulatory genomic regions encompassing 2,048,698 CpG sites. FINDINGS: We find and replicate 3 previously unreported T2D differentially methylated CpG positions (T2D-DMPs) at FDR 5% in RGL3, NGB and OTX2, and 20 signals at FDR 25%, of which 14 replicated. Integrating genetic variation and T2D-discordant monozygotic twin analyses, we identify both genetic-based and genetic-independent T2D-DMPs. The signals annotate to genes with established GWAS and EWAS links to T2D and its complications, including blood pressure (RGL3) and eye disease (OTX2). INTERPRETATION: The results help to improve our understanding of T2D disease pathogenesis and progression and may provide biomarkers for its complications. FUNDING: Funding acknowledgements for each cohort can be found in the Supplementary Note.

Conserved signatures of the canine faecal microbiome are associated with metronidazole treatment and recovery.

Antibiotic resistance is recognised as one of the biggest global threats to human and animal health. Understanding the influence of antibiotics on the canine microbiome is important to know the potential mid-to-long term effects on dysbiosis and mitigate side-effects such as antibiotic-associated diarrhoea. In this study, metronidazole was prescribed to 22 dogs for suspected giardiasis after exhibiting gastrointestinal symptoms such as diarrhoea and/or vomiting. Faecal samples were collected before, during seven days of treatment, and six months post-cessation. Faecal microbiota was assessed with 16S rRNA sequencing. Shannon diversity was reduced for up to three days after the treatment ended, and an altered community persisted for four to six weeks. All dogs recovered to a similar microbiome composition as pre-treatment. Immediately after receiving metronidazole, an increase in the relative abundance of the genera Lactobacillus, Bifidobacterium, and Enterococcus was observed. This may be due to antibiotic resistance commonly exhibited by these organisms. One-to-two weeks post-cessation, several other genera that were sensitive to the antibiotic recovered in abundances, with taxa belonging to the Erysipelotrichaceae family particularly driving composition change. Many of the bacteria initially reduced were associated with carbohydrate fermentation. This suggests scope exists to explore interventions to augment gastrointestinal health and support the re-establishment of the microbiome.

Genetic impacts on DNA methylation help elucidate regulatory genomic processes.

BACKGROUND: Pinpointing genetic impacts on DNA methylation can improve our understanding of pathways that underlie gene regulation and disease risk. RESULTS: We report heritability and methylation quantitative trait locus (meQTL) analysis at 724,499 CpGs profiled with the Illumina Infinium MethylationEPIC array in 2358 blood samples from three UK cohorts. Methylation levels at 34.2% of CpGs are affected by SNPs, and 98% of effects are cis-acting or within 1 Mbp of the tested CpG. Our results are consistent with meQTL analyses based on the former Illumina Infinium HumanMethylation450 array. Both SNPs and CpGs with meQTLs are overrepresented in enhancers, which have improved coverage on this platform compared to previous approaches. Co-localisation analyses across genetic effects on DNA methylation and 56 human traits identify 1520 co-localisations across 1325 unique CpGs and 34 phenotypes, including in disease-relevant genes, such as USP1 and DOCK7 (total cholesterol levels), and ICOSLG (inflammatory bowel disease). Enrichment analysis of meQTLs and integration with expression QTLs give insights into mechanisms underlying cis-meQTLs (e.g. through disruption of transcription factor binding sites for CTCF and SMC3) and trans-meQTLs (e.g. through regulating the expression of ACD and SENP7 which can modulate DNA methylation at distal sites). CONCLUSIONS: Our findings improve the characterisation of the mechanisms underlying DNA methylation variability and are informative for prioritisation of GWAS variants for functional follow-ups. The MeQTL EPIC Database and viewer are available online at https://epicmeqtl.kcl.ac.uk .

Integrative genomic analyses in adipocytes implicate DNA methylation in human obesity and diabetes.

DNA methylation variations are prevalent in human obesity but evidence of a causative role in disease pathogenesis is limited. Here, we combine epigenome-wide association and integrative genomics to investigate the impact of adipocyte DNA methylation variations in human obesity. We discover extensive DNA methylation changes that are robustly associated with obesity (N = 190 samples, 691 loci in subcutaneous and 173 loci in visceral adipocytes, P < 1 × 10-7). We connect obesity-associated methylation variations to transcriptomic changes at >500 target genes, and identify putative methylation-transcription factor interactions. Through Mendelian Randomisation, we infer causal effects of methylation on obesity and obesity-induced metabolic disturbances at 59 independent loci. Targeted methylation sequencing, CRISPR-activation and gene silencing in adipocytes, further identifies regional methylation variations, underlying regulatory elements and novel cellular metabolic effects. Our results indicate DNA methylation is an important determinant of human obesity and its metabolic complications, and reveal mechanisms through which altered methylation may impact adipocyte functions.

Multi-omics analyses demonstrate a critical role for EHMT1 methyltransferase in transcriptional repression during oogenesis.

EHMT1 (also known as GLP) is a multifunctional protein, best known for its role as an H3K9me1 and H3K9me2 methyltransferase through its reportedly obligatory dimerization with EHMT2 (also known as G9A). Here, we investigated the role of EHMT1 in the oocyte in comparison to EHMT2 using oocyte-specific conditional knockout mouse models (Ehmt2 cKO, Ehmt1 cKO, Ehmt1/2 cDKO), with ablation from the early phase of oocyte growth. Loss of EHMT1 in Ehmt1 cKO and Ehmt1/2 cDKO oocytes recapitulated meiotic defects observed in the Ehmt2 cKO; however, there was a significant impairment in oocyte maturation and developmental competence in Ehmt1 cKO and Ehmt1/2 cDKO oocytes beyond that observed in the Ehmt2 cKO. Consequently, loss of EHMT1 in oogenesis results, upon fertilization, in mid-gestation embryonic lethality. To identify H3K9 methylation and other meaningful biological changes in each mutant to explore the molecular functions of EHMT1 and EHMT2, we performed immunofluorescence imaging, multi-omics sequencing, and mass spectrometry (MS)-based proteome analyses in cKO oocytes. Although H3K9me1 was depleted only upon loss of EHMT1, H3K9me2 was decreased, and H3K9me2-enriched domains were eliminated equally upon loss of EHMT1 or EHMT2. Furthermore, there were more significant changes in the transcriptome, DNA methylome, and proteome in Ehmt1/2 cDKO than Ehmt2 cKO oocytes, with transcriptional derepression leading to increased protein abundance and local changes in genic DNA methylation in Ehmt1/2 cDKO oocytes. Together, our findings suggest that EHMT1 contributes to local transcriptional repression in the oocyte, partially independent of EHMT2, and is critical for oogenesis and oocyte developmental competence.

Adipose methylome integrative-omic analyses reveal genetic and dietary metabolic health drivers and insulin resistance classifiers.

BACKGROUND: There is considerable evidence for the importance of the DNA methylome in metabolic health, for example, a robust methylation signature has been associated with body mass index (BMI). However, visceral fat (VF) mass accumulation is a greater risk factor for metabolic disease than BMI alone. In this study, we dissect the subcutaneous adipose tissue (SAT) methylome signature relevant to metabolic health by focusing on VF as the major risk factor of metabolic disease. We integrate results with genetic, blood methylation, SAT gene expression, blood metabolomic, dietary intake and metabolic phenotype data to assess and quantify genetic and environmental drivers of the identified signals, as well as their potential functional roles. METHODS: Epigenome-wide association analyses were carried out to determine visceral fat mass-associated differentially methylated positions (VF-DMPs) in SAT samples from 538 TwinsUK participants. Validation and replication were performed in 333 individuals from 3 independent cohorts. To assess functional impacts of the VF-DMPs, the association between VF and gene expression was determined at the genes annotated to the VF-DMPs and an association analysis was carried out to determine whether methylation at the VF-DMPs is associated with gene expression. Further epigenetic analyses were carried out to compare methylation levels at the VF-DMPs as the response variables and a range of different metabolic health phenotypes including android:gynoid fat ratio (AGR), lipids, blood metabolomic profiles, insulin resistance, T2D and dietary intake variables. The results from all analyses were integrated to identify signals that exhibit altered SAT function and have strong relevance to metabolic health. RESULTS: We identified 1181 CpG positions in 788 genes to be differentially methylated with VF (VF-DMPs) with significant enrichment in the insulin signalling pathway. Follow-up cross-omic analysis of VF-DMPs integrating genetics, gene expression, metabolomics, diet, and metabolic traits highlighted VF-DMPs located in 9 genes with strong relevance to metabolic disease mechanisms, with replication of signals in FASN, SREBF1, TAGLN2, PC and CFAP410. PC methylation showed evidence for mediating effects of diet on VF. FASN DNA methylation exhibited putative causal effects on VF that were also strongly associated with insulin resistance and methylation levels in FASN better classified insulin resistance (AUC=0.91) than BMI or VF alone. CONCLUSIONS: Our findings help characterise the adiposity-associated methylation signature of SAT, with insights for metabolic disease risk.

Pulmonary Function and Blood DNA Methylation: A Multiancestry Epigenome-Wide Association Meta-analysis.

Rationale: Methylation integrates factors present at birth and modifiable across the lifespan that can influence pulmonary function. Studies are limited in scope and replication. Objectives: To conduct large-scale epigenome-wide meta-analyses of blood DNA methylation and pulmonary function. Methods: Twelve cohorts analyzed associations of methylation at cytosine-phosphate-guanine probes (CpGs), using Illumina 450K or EPIC/850K arrays, with FEV1, FVC, and FEV1/FVC. We performed multiancestry epigenome-wide meta-analyses (total of 17,503 individuals; 14,761 European, 2,549 African, and 193 Hispanic/Latino ancestries) and interpreted results using integrative epigenomics. Measurements and Main Results: We identified 1,267 CpGs (1,042 genes) differentially methylated (false discovery rate, <0.025) in relation to FEV1, FVC, or FEV1/FVC, including 1,240 novel and 73 also related to chronic obstructive pulmonary disease (1,787 cases). We found 294 CpGs unique to European or African ancestry and 395 CpGs unique to never or ever smokers. The majority of significant CpGs correlated with nearby gene expression in blood. Findings were enriched in key regulatory elements for gene function, including accessible chromatin elements, in both blood and lung. Sixty-nine implicated genes are targets of investigational or approved drugs. One example novel gene highlighted by integrative epigenomic and druggable target analysis is TNFRSF4. Mendelian randomization and colocalization analyses suggest that epigenome-wide association study signals capture causal regulatory genomic loci. Conclusions: We identified numerous novel loci differentially methylated in relation to pulmonary function; few were detected in large genome-wide association studies. Integrative analyses highlight functional relevance and potential therapeutic targets. This comprehensive discovery of potentially modifiable, novel lung function loci expands knowledge gained from genetic studies, providing insights into lung pathogenesis.

Epigenetic changes induced by in utero dietary challenge result in phenotypic variability in successive generations of mice.

Transmission of epigenetic information between generations occurs in nematodes, flies and plants, mediated by specialised small RNA pathways, modified histones and DNA methylation. Similar processes in mammals can also affect phenotype through intergenerational or trans-generational mechanisms. Here we generate a luciferase knock-in reporter mouse for the imprinted Dlk1 locus to visualise and track epigenetic fidelity across generations. Exposure to high-fat diet in pregnancy provokes sustained re-expression of the normally silent maternal Dlk1 in offspring (loss of imprinting) and increased DNA methylation at the somatic differentially methylated region (sDMR). In the next generation heterogeneous Dlk1 mis-expression is seen exclusively among animals born to F1-exposed females. Oocytes from these females show altered gene and microRNA expression without changes in DNA methylation, and correct imprinting is restored in subsequent generations. Our results illustrate how diet impacts the foetal epigenome, disturbing canonical and non-canonical imprinting mechanisms to modulate the properties of successive generations of offspring.

Leptin Signaling in the Ovary of Diet-Induced Obese Mice Regulates Activation of NOD-Like Receptor Protein 3 Inflammasome.

Obesity leads to ovarian dysfunction and the establishment of local leptin resistance. The aim of our study was to characterize the levels of NOD-like receptor protein 3 (NLRP3) inflammasome activation in ovaries and liver of mice during obesity progression. Furthermore, we tested the putative role of leptin on NLRP3 regulation in those organs. C57BL/6J female mice were treated with equine chorionic gonadotropin (eCG) or human chorionic gonadotropin (hCG) for estrous cycle synchronization and ovary collection. In diet-induced obesity (DIO) protocol, mice were fed chow diet (CD) or high-fat diet (HFD) for 4 or 16 weeks, whereas in the hyperleptinemic model (LEPT), mice were injected with leptin for 16 days (16 L) or saline (16 C). Finally, the genetic obese leptin-deficient ob/ob (+/? and -/-) mice were fed CD for 4 week. Either ovaries and liver were collected, as well as cumulus cells (CCs) after superovulation from DIO and LEPT. The estrus cycle synchronization protocol showed increased protein levels of NLRP3 and interleukin (IL)-18 in diestrus, with this stage used for further sample collections. In DIO, protein expression of NLRP3 inflammasome components was increased in 4 week HFD, but decreased in 16 week HFD. Moreover, NLRP3 and IL-1ß were upregulated in 16 L and downregulated in ob/ob. Transcriptome analysis of CC showed common genes between LEPT and 4 week HFD modulating NLRP3 inflammasome. Liver analysis showed NLRP3 protein upregulation after 16 week HFD in DIO, but also its downregulation in ob/ob-/-. We showed the link between leptin signaling and NLRP3 inflammasome activation in the ovary throughout obesity progression in mice, elucidating the molecular mechanisms underpinning ovarian failure in maternal obesity.

Genomic and phenotypic insights from an atlas of genetic effects on DNA methylation.

Characterizing genetic influences on DNA methylation (DNAm) provides an opportunity to understand mechanisms underpinning gene regulation and disease. In the present study, we describe results of DNAm quantitative trait locus (mQTL) analyses on 32,851 participants, identifying genetic variants associated with DNAm at 420,509 DNAm sites in blood. We present a database of >270,000 independent mQTLs, of which 8.5% comprise long-range (trans) associations. Identified mQTL associations explain 15-17% of the additive genetic variance of DNAm. We show that the genetic architecture of DNAm levels is highly polygenic. Using shared genetic control between distal DNAm sites, we constructed networks, identifying 405 discrete genomic communities enriched for genomic annotations and complex traits. Shared genetic variants are associated with both DNAm levels and complex diseases, but only in a minority of cases do these associations reflect causal relationships from DNAm to trait or vice versa, indicating a more complex genotype-phenotype map than previously anticipated.

Identical twins carry a persistent epigenetic signature of early genome programming.

Monozygotic (MZ) twins and higher-order multiples arise when a zygote splits during pre-implantation stages of development. The mechanisms underpinning this event have remained a mystery. Because MZ twinning rarely runs in families, the leading hypothesis is that it occurs at random. Here, we show that MZ twinning is strongly associated with a stable DNA methylation signature in adult somatic tissues. This signature spans regions near telomeres and centromeres, Polycomb-repressed regions and heterochromatin, genes involved in cell-adhesion, WNT signaling, cell fate, and putative human metastable epialleles. Our study also demonstrates a never-anticipated corollary: because identical twins keep a lifelong molecular signature, we can retrospectively diagnose if a person was conceived as monozygotic twin.

Estrogen and COVID-19 symptoms: Associations in women from the COVID Symptom Study.

It has been widely observed that adult men of all ages are at higher risk of developing serious complications from COVID-19 when compared with women. This study aimed to investigate the association of COVID-19 positivity and severity with estrogen exposure in women, in a population based matched cohort study of female users of the COVID Symptom Study application in the UK. Analyses included 152,637 women for menopausal status, 295,689 women for exogenous estrogen intake in the form of the combined oral contraceptive pill (COCP), and 151,193 menopausal women for hormone replacement therapy (HRT). Data were collected using the COVID Symptom Study in May-June 2020. Analyses investigated associations between predicted or tested COVID-19 status and menopausal status, COCP use, and HRT use, adjusting for age, smoking and BMI, with follow-up age sensitivity analysis, and validation in a subset of participants from the TwinsUK cohort. Menopausal women had higher rates of predicted COVID-19 (P = 0.003). COCP-users had lower rates of predicted COVID-19 (P = 8.03E-05), with reduction in hospital attendance (P = 0.023). Menopausal women using HRT or hormonal therapies did not exhibit consistent associations, including increased rates of predicted COVID-19 (P = 2.22E-05) for HRT users alone. The findings support a protective effect of estrogen exposure on COVID-19, based on positive association between predicted COVID-19 with menopausal status, and negative association with COCP use. HRT use was positively associated with COVID-19, but the results should be considered with caution due to lack of data on HRT type, route of administration, duration of treatment, and potential unaccounted for confounders and comorbidities.

Childhood growth and development and DNA methylation age in mid-life.

BACKGROUND: In the first study of its kind, we examine the association between growth and development in early life and DNAm age biomarkers in mid-life. METHODS: Participants were from the Medical Research Council National Survey of Health and Development (n = 1376). Four DNAm age acceleration (AgeAccel) biomarkers were measured when participants were aged 53 years: AgeAccelHannum; AgeAccelHorvath; AgeAccelLevine; and AgeAccelGrim. Exposure variables included: relative weight gain (standardised residuals from models of current weight z-score on current height, and previous weight and height z-scores); and linear growth (standardised residuals from models of current height z-score on previous height and weight z-scores) during infancy (0-2 years, weight gain only), early childhood (2-4 years), middle childhood (4-7 years) and late childhood to adolescence (7-15 years); age at menarche; and pubertal stage for men at 14-15 years. The relationship between relative weight gain and linear growth and AgeAccel was investigated using conditional growth models. We replicated analyses from the late childhood to adolescence period and pubertal timing among 240 participants from The National Child and Development Study (NCDS). RESULTS: A 1SD increase in relative weight gain in late childhood to adolescence was associated with 0.50 years (95% CI 0.20, 0.79) higher AgeAccelGrim. Although the CI includes the null, the estimate was similar in NCDS [0.57 years (95% CI - 0.01, 1.16)] There was no strong evidence that relative weight gain and linear growth in childhood was associated with any other AgeAccel biomarker. There was no relationship between pubertal timing in men and AgeAccel biomarkers. Women who reached menarche ≥ 12 years had 1.20 years (95% CI 0.15, 2.24) higher AgeAccelGrim on average than women who reached menarche < 12 years; however, this was not replicated in NCDS and was not statistically significant after Bonferroni correction. CONCLUSIONS: Our findings generally do not support an association between growth and AgeAccel biomarkers in mid-life. However, we found rapid weight gain during pubertal development, previously related to higher cardiovascular disease risk, to be associated with older AgeAccelGrim. Given this is an exploratory study, this finding requires replication.

Epigenome-wide association meta-analysis of DNA methylation with coffee and tea consumption.

Coffee and tea are extensively consumed beverages worldwide which have received considerable attention regarding health. Intake of these beverages is consistently linked to, among others, reduced risk of diabetes and liver diseases; however, the mechanisms of action remain elusive. Epigenetics is suggested as a mechanism mediating the effects of dietary and lifestyle factors on disease onset. Here we report the results from epigenome-wide association studies (EWAS) on coffee and tea consumption in 15,789 participants of European and African-American ancestries from 15 cohorts. EWAS meta-analysis of coffee consumption reveals 11 CpGs surpassing the epigenome-wide significance threshold (P-value <1.1×10-7), which annotated to the AHRR, F2RL3, FLJ43663, HDAC4, GFI1 and PHGDH genes. Among them, cg14476101 is significantly associated with expression of the PHGDH and risk of fatty liver disease. Knockdown of PHGDH expression in liver cells shows a correlation with expression levels of genes associated with circulating lipids, suggesting a role of PHGDH in hepatic-lipid metabolism. EWAS meta-analysis on tea consumption reveals no significant association, only two CpGs annotated to CACNA1A and PRDM16 genes show suggestive association (P-value <5.0×10-6). These findings indicate that coffee-associated changes in DNA methylation levels may explain the mechanism of action of coffee consumption in conferring risk of diseases.

Revisiting the Impact of Local Leptin Signaling in Folliculogenesis and Oocyte Maturation in Obese Mothers.

The complex nature of folliculogenesis regulation accounts for its susceptibility to maternal physiological fitness. In obese mothers, progressive expansion of adipose tissue culminates with severe hyperestrogenism and hyperleptinemia with detrimental effects for ovarian performance. Indeed, maternal obesity is associated with the establishment of ovarian leptin resistance. This review summarizes current knowledge on potential effects of impaired leptin signaling throughout folliculogenesis and oocyte developmental competence in mice and women.

Life course socioeconomic position and DNA methylation age acceleration in mid-life.

BACKGROUND: Ageing biomarkers can help us better understand how well-established socioeconomic position (SEP) disparities in ageing occur. A promising new set of DNAm methylation (DNAm)-based ageing biomarkers indicate through their age acceleration (AA) measures if biological ageing is slower or faster than chronological ageing. Few studies have investigated the association between SEP and DNAm AA. METHODS: We used linear regression to examine the sex-adjusted relationships between childhood social class, adult social class, intergenerational social class change, education and adult household earnings with first (Horvath AA and Hannum AA) and second generation (PhenoAge AA and GrimAge AA) DNAm AA markers using data from the MRC National Survey of Health and Development. RESULTS: In the first-generation biomarkers, there was little evidence of any associations with Horvath AA but associations of childhood social class and income with Hannum AA were observed. Strong associations were seen between greater disadvantage in childhood and adult SEP and greater AA in the second generation biomarkers. For example, those with fathers in an unskilled occupational social class in childhood had 3.6 years greater PhenoAge AA (95% CI 1.8 to 5.4) than those with fathers from a professional social class. Individuals without qualifications had higher AA compared with those with higher education (4.1 years greater GrimAge AA (95% CI 3.1 to 5.0)). CONCLUSION: Our findings highlight the importance of exposure to social disadvantage in childhood to the biological ageing process. The second generation clocks appear to be more sensitive to the accumulation of social disadvantage across the life course.

Epigenome-Wide Association Study of Thyroid Function Traits Identifies Novel Associations of fT3 With KLF9 and DOT1L.

CONTEXT: Circulating concentrations of free triiodothyronine (fT3), free thyroxine (fT4), and thyrotropin (TSH) are partly heritable traits. Recent studies have advanced knowledge of their genetic architecture. Epigenetic modifications, such as DNA methylation (DNAm), may be important in pituitary-thyroid axis regulation and action, but data are limited. OBJECTIVE: To identify novel associations between fT3, fT4, and TSH and differentially methylated positions (DMPs) in the genome in subjects from 2 Australian cohorts. METHOD: We performed an epigenome-wide association study (EWAS) of thyroid function parameters and DNAm using participants from: Brisbane Systems Genetics Study (median age 14.2 years, n = 563) and the Raine Study (median age 17.0 years, n = 863). Plasma fT3, fT4, and TSH were measured by immunoassay. DNAm levels in blood were assessed using Illumina HumanMethylation450 BeadChip arrays. Analyses employed generalized linear mixed models to test association between DNAm and thyroid function parameters. Data from the 2 cohorts were meta-analyzed. RESULTS: We identified 2 DMPs with epigenome-wide significant (P < 2.4E-7) associations with TSH and 6 with fT3, including cg00049440 in KLF9 (P = 2.88E-10) and cg04173586 in DOT1L (P = 2.09E-16), both genes known to be induced by fT3. All DMPs had a positive association between DNAm and TSH and a negative association between DNAm and fT3. There were no DMPs significantly associated with fT4. We identified 23 differentially methylated regions associated with fT3, fT4, or TSH. CONCLUSIONS: This study has demonstrated associations between blood-based DNAm and both fT3 and TSH. This may provide insight into mechanisms underlying thyroid hormone action and/or pituitary-thyroid axis function.

Identification of a unique epigenetic profile in women with diminished ovarian reserve.

OBJECTIVE: To investigate whether epigenetic profiles of mural granulosa cells (MGC) and leukocytes from women with diminished ovarian reserve (DOR) differ from those of women with normal or high ovarian reserve. DESIGN: Prospectively collected material from a multicenter cohort of women undergoing fertility treatment. SETTING: Private and university-based facilities for clinical services and research. PATIENT(S): One hundred and nineteen women of various ages and ovarian reserve status (antimüllerian hormone level) who provided blood samples and MGC. INTERVENTION(S): None. MAIN OUTCOME MEASURE(S): Measures of epigenetic aging rates from whole-genome methylation array data: DNA methylation variability, age acceleration, DNA methylation telomere length estimator (DNAmTL), and accumulation of epimutations. RESULT(S): Comparison of DOR or high ovarian reserve samples to controls (normal ovarian reserve) showed differential methylation variability between DOR and normal samples at 4,199 CpGs in MGC, and 447 between high and normal (false-discovery rate < 0.05). Variable sites in MGC from DOR were enriched in regions marked with the repressive histone modification H3K27me3, and also included genes involved in folliculogenesis, such as insulin growth factor 2 (IGF2) and antimüllerian hormone (AMH). Regardless of ovarian reserve, very few signals were detected in leukocytes, and no overlaps with those in MGC were found. Furthermore, we found a higher number of epimutations in MGC from women with DOR (Kruskal-Wallis test, difference in mean = 3,485). CONCLUSION(S): The somatic cells of human ovarian follicles have a distinctive epigenetic profile in women with DOR. A high frequency of epimutations suggests premature aging. Ovarian reserve status was not reflected in the leukocyte epigenetic profile.

Increased transcriptome variation and localised DNA methylation changes in oocytes from aged mice revealed by parallel single-cell analysis.

Advancing maternal age causes a progressive reduction in fertility. The decline in developmental competence of the oocyte with age is likely to be a consequence of multiple contributory factors. Loss of epigenetic quality of the oocyte could impair early developmental events or programme adverse outcomes in offspring that manifest only later in life. Here, we undertake joint profiling of the transcriptome and DNA methylome of individual oocytes from reproductively young and old mice undergoing natural ovulation. We find reduced complexity as well as increased variance in the transcriptome of oocytes from aged females. This transcriptome heterogeneity is reflected in the identification of discrete sub-populations. Oocytes with a transcriptome characteristic of immature chromatin configuration (NSN) clustered into two groups: one with reduced developmental competence, as indicated by lower expression of maternal effect genes, and one with a young-like transcriptome. Oocytes from older females had on average reduced CpG methylation, but the characteristic bimodal methylation landscape of the oocyte was preserved. Germline differentially methylated regions of imprinted genes were appropriately methylated irrespective of age. For the majority of differentially expressed transcripts, the absence of correlated methylation changes suggests a post-transcriptional basis for most age-related effects on the transcriptome. However, we did find differences in gene body methylation at which there were corresponding changes in gene expression, indicating age-related effects on transcription that translate into methylation differences. Interestingly, oocytes varied in expression and methylation of these genes, which could contribute to variable competence of oocytes or penetrance of maternal age-related phenotypes in offspring.

Landscape of Genome-Wide DNA Methylation of Colorectal Cancer Metastasis.

Colorectal cancer is a heterogeneous disease caused by both genetic and epigenetics factors. Analysing DNA methylation changes occurring during colorectal cancer progression and metastasis formation is crucial for the identification of novel epigenetic markers of patient prognosis. Genome-wide methylation sequencing of paired samples of colon (normal adjacent, primary tumour and lymph node metastasis) showed global hypomethylation and CpG island (CGI) hypermethylation of primary tumours compared to normal. In metastasis we observed high global and non-CGI regions methylation, but lower CGI methylation, compared to primary tumours. Gene ontology analysis showed shared biological processes between hypermethylated CGIs in metastasis and primary tumours. After complementary analysis with The Cancer Genome Atlas (TCGA) cohort, FIGN, HTRA3, BDNF, HCN4 and STAC2 genes were found associated with poor survival. We mapped the methylation landscape of colon normal tissues, primary tumours and lymph node metastasis, being capable of identified methylation changes throughout the genome. Furthermore, we found five genes with potential for methylation biomarkers of poor prognosis in colorectal cancer patients.